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Abstract. The phononlike excitations of (anti-)instanton (II) liquid due 
to adiabatic variations of vacuum wave functions are studied in this pa- 
per. The kinetic energy term is found and the proper effective Lagrangian 
for such excitations is evaluated. The properties of their spectrum, while 
corresponding masses are defined by Aqcd, are investigated. 



The model of (anti-)instanton liquid correctly seizes many nonpertur- 
bative phenomena and important vacuum features such as chiral symmetry 
breaking, the presence of gluon condensate and topological susceptibility 
[1],[2]. It is usually supposed that the corresponding functional integral in 
this approach is saturated by quasi-classical configurations close to the ex- 
act solutions of the Yang- Mills equations (the Euclidean solutions called the 
(anti-)instantons) and the wave function of vacuum, being homogeneous in 
metric space, is properly reproduced by averaging over their collective coor- 
dinates. In the II liquid approach one takes the superposition ansatz of the 
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pseudo-particle (PP) fields as one of the simplest relevant approximations 
to the 'genuine' vacuum configuration 

N 

A„{x) =J2M x ^i) • (1) 
i=i 

Here A^(x;~fi) denotes the field of a singled (anti-)instanton in singular 
gauge with 4N C (for the SU(N C ) group) coordinates 7 = (p, z, fl), of size 
p with the coordinate of its centre z, Q as its colour orientation and 

A«(x; 1 ) = -n a % tMU ^^-^ , y = x-z, (2) 

where 77 is the 't Hooft symbol [3]. For anti-instanton 77 — > 77 (making the 
choice of the singular gauge allows us to sum up the solution preserving 
the asymptotic behaviour). For simplicity we shall not introduce different 
symbols for instanton and anti-instanton, and then in the superposition of 
Eq.(l) A?" implies the PP total number in the 4-volume V system with the 
density n = N/V . The action of the instanton liquid model is introduced 
by the following functional 

(S) = Jd A zJ dp n{p)s{p) . (3) 

The integration should be performed over the system volume along with 
averaging the one instanton action s(p) weighted by instanton size distribu- 
tion function n(p). Then an action per one instanton is given by well-known 
expression 

Sl (p) = P(p) + 51n(Ap) - ln£ 2JV <= + 0f l p 2 J d Pl n(pi)p\ , (4) 

with the Gell-Mann-Low beta function f3(p) = — lnC7v c — 61n(Ap) , A = 
Ajfg = 0.92Ap.y. , and constant Cn c depending on the regularization 

scheme, here Cn c ~ t-^(AT C ^i)!(iV 6 ^2)! ' ^ = T^ c ' anc ^ P arame t ers 
(5 = 0{p) and (3 = (5 + In Cn c are the (3 function values at the fixed quantity 
of p (average instanton size). 

Some terms (inperfect p dependence) of Eq.(4) could be obtained in 
the classical Yang-Mills theory with one-loop (quantum) corrections taken 
into account and resulting in a modification of coupling constant g at the 
distinct scales. Indeed, the first term is the one instanton action 8tt 2 /g 2 
with the p-dependence of g corrected. The last term of Eq.(4) describes 
the pair interaction of PP in the instanton ensemble with the constant 
£ characterizing, in a sense, the intensity of interaction £ 2 = xiv^T 71 " 2, 
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The smallness of characteristic instanton liquid parameter np 4 ('packing 
fraction') allows us to drop out the p dependence of the /3-function. The 
logarithmic terms correspond to the pre-exponential factor contribution to 
the functional integral and are of a genuine quantum nature. In the second 
term the p 5 factor makes the integration elements dp and d 4 z dimensionless. 
Finally, the third term is logarithm of a square root of the one instanton 
action f3 raised to the power AN C . The latter is just the zero mode number 
of the one instanton solution and the corresponding p dependence may be 
again omitted because of small logarithmic contribution. 

Taking the exponential form for the distribution function over the in- 
stanton action n(p) ~ e~ Sl ^ p \ we obtain directly from Eq.(3) the self- 
consistent description of equilibrium state of instanton liquid with the well- 
known ground state p(p) = p- 5 p 2N <=e~ ,3 ^~ up ' 2 '/ p2 , where 
u=(b-4)/2, (^) 2 = ^, n = jdpn(p), J = J dp p>n{p)/n. 

This argument corresponds to the maximum principle of [2]. Approach- 
ing the functional of the averaged action per one instanton as a local form 
( Sl ) = J dp s 1 (p)n(p)/n where si(p) = /3(p) + 51n(Ap)-ln/3 27Vc + /3^ 2 p 2 np2, 
using n(p) = Ce~ s( - pS) (C = const) as a distribution function (actually it 
makes the problem self-consistent because an equilibrium distribution func- 
tion should be dependent on an action only) and taking the variation of 
(s) — (si) = J dp {s(p) — si(p)}e~ s ( p ^ jn over s(p) one may come to the 
result s(p) = si(p) + const keeping into the mind an arbitrary normaliza- 
tion. Then maximizing the action per one instanton over p, for example, we 
obtain an explicit description (E.D.) of ground state of the instanton liquid 

through the equilibrium value of p = e 2v ~ 1 , close to the values obtained 
in [2] (D.P.) 
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where Nf is the number of flavours, = 3. 

The distribution p(p) has obvious physical meaning, namely, the quan- 
tity d 4 z dp p(p) is proportional to the probability to find an instanton of 
size p at some point of a volume element d 4 z. At small p the behaviour of 
distribution function is stipulated by the quantum-mechanical uncertainty 
principle preventing a solution being compressed at a point (radiative cor- 
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rection). At large p the constraint comes from the repulsive interaction 
between the PP which is amplified with (anti-)instanton size increasing. 

Deriving Eq.(3) we should average over the instanton positions in a 
metric space. It is clear that the characteristic size, which has to be taken 
into account, should be larger enough than the mean instanton size p. But 
it should not be too large because the far ranged fragments of instanton 
liquid are not 'causally' dependent. The vacuum wave function is expected 
to be homogeneous on this scale L > R (R is an average distance between 
the PP). Let us remind that each PP contributes to the functional integral 
with the weight factor proportional to 1/V, V = L 4 . The characteristic 
configuration saturating the functional integral is taken as the superposition 
(1) with N pseudoparticles in the volume V. If one supposes that the PP 
number in an ensemble is still appropriate to consider them separately, then 
denoting AN(pi) as the PP number of size p € (pi, pi + Ap), K as the 
number of partitions within the interval (pi, pf), Eq.(l) may be rewritten 
in the following form 

K AN(pi) 

M X )=Y, E M x Mj) > ( 5 ) 

i=l j=l 

where A^(x;i,jj) is the (anti-)instanton solution with the calibrated size 
and 7 = (z, f2) stands for the coordinate of its centre and colour orienta- 
tion. By definition Y^f=i AN(pi) = N. Further, introducing the distribution 
function n(p) = A ^ p ^ y , and normalizing it as Ya=i n (Pi) A/? = A^ (in 
the continual limit Ap — ^ it is valid V J dp n(p) = N) one can calculate 
the classical action S c = \ J d A x G 2 ^ v of this configuration averaging over 
the instanton positions in the metric and colour spaces. As a result (with 
the superposition ansatz Eq.(l)) one instanton actions and the PP pair 
interactions only contribute to the average system action 

( S c) = fl / ^r dn * S c=J d 4 z I dpn(p) j^f + / d Pl n{Pi)p\ 

where dVL is a measure in the colour space with the unit normalization. As 
above mentioned (S c ) including one loop corrections will then contribute 
to the functional integral. 

It is easy to understand that Eq.(3) describes properly even non-equi- 
librium states of the instanton liquid when the distribution function n(p) 
does not coincide with the ground state one p(p). Moreover, it allows us 
to generalize Eq.(3) for the non- homogeneous liquid, when the size of the 
non-homogeneity obeys the obvious constraint A > L > p. 

In what follows, we study the excitations of II liquid induced by adi- 
abatic dilatational deformations of the instanton solutions. Then, as the 
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configurations saturating the functional integral we consider not the in- 
stanton solution itself but the quasizero modes which are parametrically 
very close in the functional space (a direction does exist where the action 
varies slowly) to the zero modes. The guiding idea of selecting a deforma- 
tion originates from transparent observation. The deformations measured 
in units of the action (here q, p are the generalized coordinate and 

momentum) have a physical meaning only. However, the instantons are 
characterized by 'static' coordinates 7 and, therefore, need to appoint the 
conjugated momenta. It looks quite natural for the variable, for example, 
p to introduce those as p = dp/dx^. 

Let us calculate first of all the corrections for the one-instanton action. 
Dealing with superposition ansatz Eq.(l) again, one should include addi- 
tional contribution to the chromoelectric field 

G% = G%, + <fa,. (6) 

with the first term of strength tensor (s.t.) corresponding to the contribu- 
tion generated by the instanton profile 

r* 8 p2 (-- 1 — ypVp_- IMh 
g {y 2 + p 2 ) 2 \2 natxv avp y 2 T]aixp y 2 

and in adiabatic approximation the corrections have the form 

QA a 4 y v p 

dti = ~q^P = ~Vaiv ^ 2 + ^2)2 P-< 9?j = > 9i4 = -gti , % i3 = 1) 2, 3 . 

The adiabatic constraint g^ u <C G^ means that the variation of instanton 
size is much smaller than the characteristic transformation scale of the PP 
field, p <C 0(1). Then calculating the corrections for the action, it is reason- 
able to take out p beyond the integral and the one instanton contribution 
to the action turns out to be 



d x G^ ~ — + C p + — p , (7) 



where p should be taken as the mean rate of slow solution deformation at a 
characteristic instanton lifetime ~ p. For simplicity, one may take it in the 
centre of the instanton p(z) = p(x, z)\ x=z . The constant C = (because the 
first variation of the action SS/5A for the solution itself equals to zero). For 
the 'kinematicaP K-term we have K s ,t. = ^jp-- The overt p dependence of k 
is lacking because of the scale invariance. It arises with the renormalization 
of the coupling constant (in a regular gauge the result is the same) . Being 
within the ansatz (1) we have considered only the corrections induced by 
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the variation of strength tensors, but not those resulting from a possible 
variation of fields (2). Bearing in mind the form of potentials in regular 
(r.g.) and singular (s.g.) gauges 

= \»„ d^Hy 2 + p 2 ) ( r.g.) ; A® = ~fj a ^ d„ln ^1 + ( s.g.) 

we find the adiabatic corrections = + as follows: 

2 p 

= -f] ail A 2 , 2 P (r.g.) • (8) 

^ g y 2 + p 2 

The substitution r] — > — fj brings about the transition from regular gauge 
to a singular one. Using the admixture of chromoelectric and chromo- 
magnetic fields generated by the to saturate the functional integral as 
in Eq.(6) we drop the terms of higher order than O(p) out. Thus, we come 
to the result for the 'kinematicaP term re as 1 

32vr 2 . , 
« = -jr ( r -S-) • ( 9 ) 

But within the accuracy taken at calculating Eq.(3) we are permitted to 
fix re at some point as re(p). 

Analysing the variations of the functional integration result when having 
calculated for the quasizero mode in the adiabatic approximation we see 
that because of the kinetic energy smallness it is certainly permitted to 
neglect its impact on the pre-exponential factor, which is small itself. The 
inverse influence of the pre-exponential factor on the 'kinematical' term 
is negligible as well. Thus, going the way which we have already passed 
through while calculating Eq.(3) we receive within the required order of 
accuracy 2 

(S) = Jd A zJ dp n{p) [\*P 2 + s(p)} • (10) 

It is important to remark here that the contribution of averaged PP pair in- 
teraction under the adiabatic conditions is estimated as a term with contact 
interaction 

d 4 z8vr 2 i 2 
9 2 



(S(12)) * J £ P \{z) p 2 (z) . (11) 



1 In the singular gauge we have k = | 32 % . 

2 For the case at hand the types of deformations (quasizero modes) may be simply 
counted in terms of the instanton solution. They are induced by the variations of the 
instanton size, the changes of its position and colour orientation. In fact, there is one 
more type of the quasizero modes related to two far distant instantons. 
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In particular, if the PP sizes do not vary we reproduce the well-known result 
[2] (S int (12)) = ^p- p\p\ for such a contribution. 

In the Minkowski space the factor in the curie brackets of Eq.(lO) might 
be interpreted as a mechanical system with Lagrangian C = ^ k p 2 — U e f f (p) 
if as a characteristic 'velocity' dp(x, z)jdx±\ x=z we take a deformation 
of the field p(x, z)\ x=z = p(z), dp/dx^ ~ dp/dz^, and an action per 
one instanton might be taken as a 'potential energy' U e ff(p) = (3(p) + 

5 ln(Ap) — In (3 2Nc + v £= . In the local vicinity of the potential minimum 

Pc = ^ET P 2 ' {j~dp~ = ^) * ne s yst em is oscillating and we have for the 
frequency (using the configuration corresponding to Eq.(9)) 

m 2 = — = = — = (12) 

up 2 P p 2 

while calculating the second derivative - U j f (^ L = =. 

We have only analysed the deformations in the temporary direction. 
Those in the spatial directions could be estimated by drawing the same 
arguments. Thus, the expression for the K-term keeps the form obtained 
above with the only change of rates for the appropriate gradients of function 
p(x,z), i.e. the substitution p(t) — > dp<yX '^ x ^ z ~ ^ should be performed 
for such a 'crumpled' instanton. Then the frequency of proper fluctuations 
might be interpreted as the mass term and the excitations occur to have a 
phonon-like nature 

C = \*[p 2 - VpVp] - U eff (p) , (13) 

(the cross-terms ~ / bp equal to zero identically) 3 . With the instanton liq- 
uid parameters obtained above we have for the mass term m w 1.21A. Then 
the wave length A4 ~ mT x in the X4-direction is A4A « 0.83 > V l ^K w 
0.81 > pA (the size of the phonon localization in the spatial directions can 
be arbitrary and may noticeably exceed A4, besides the number N of the 
PP forming the excitations might be pretty large). 

These numerical values obtained should be taken rather qualitatively il- 
lustrating the principle possibility to have the particle-like excitations orig- 
inated by the quasizero modes. It is clear that striving to go beyond the 
superposition ansatz one should to take into account at least a medium 
change of the instanton profile and to develop more realistic description of 

3 It is interesting to remark that the centre of instanton solution may not be shifted 
since the relevant deformations lead to the singular k (unlike the dilatational mode), 
whereas the colour coordinate variation Q gives the trivial result k = 0. 
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the instanton interactions. Certainly, what presented here essentially ex- 
ceeds the corresponding results which one could expect in the 'complete 
theory'. 

In conclusion let us emphasize that we have considered the excitations 
of the instanton liquid generated by the dilatational instanton deforma- 
tions and the adiabatic assumption leads, in principle, to a fully consistent 
picture. The quantum numbers of the excitations are the scalar 'glueball' 
one and model itself regulates the most suitable regime of such phonon- 
like deformations resulting in mass gap generation fixed by Aqcd- This 
value of the mass scale looks very reasonable and absolutely self-consistent 
for an adiabatic approximation within the instanton picture of quenched 
QCD vacuum. Indeed, it is qualitatively clear that 'hard' glueball mass 
~ 1.5 GeV (which is just known from existing theoretical investigations 
in lattice QCD and sum rule methods) certainly leads beyond our imposed 
constraints. It means that the strength of a deformation field becomes quite 
compatible with the instanton field and suppresses the PP. But the 'hard' 
glueball scale is extracted from the exponential fall of gluon correlation 
functions, and we believe that there is not well grounded to attribute this 
measurements a dynamical meaning of a signal propagation. 

Apparently, including the quark condensate an intriguing guess is to 
associate some light hadrons with these phonon-like excitations discovered 
since the preliminary evaluations of their mass spectrum look quite encour- 
aging. We believe that we could imagine them as observable in mixing with 
the quark condensat, for example, with 'sigma' meson. Moreover, the con- 
cept of the confining potential for the light quarks in the context of our 
approach seems simply irrelevant because of a stable phonon nature. 
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